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Method and device for separating acoustic signals 

The present invention relates to a method and a 
device for separating acoustic signals. 

The invention relates to the field of digital signal 
5 processing as a means of separating different acoustic 
signals from different spatial directions which are 
stereophonically picked up by two microphones at a known 
distance . 

10 The field of source separation, also referred to as 

"beam forming" is gaining in importance due to the 
increase in mobile communication as well as automatic 
processing of human speech. In very many applications, 
one problem which arises is the fact that the desired 

15 speech signal (wanted signal) is detrimentally affected 
by various types of interference. Primary examples of 
this is interference caused by background noise, 
interference from other speakers and interference from 
loudspeaker emissions of music or speech. The various 

20 types of interference require different treatments, 
depending on their nature and depending on what is known 
about the wanted signal beforehand. 

Examples of applications to ' which the invention 
25- lends itself, therefore, are communication systems in 
which the position of a speaker is known and in which 
interference occurs due to background noise or other 
speakers and loudspeaker emissions. Examples of 
applications are automotive hands-free units, in which 
30 the microphones are mounted in the rear-view mirror, for 
example, and a so-called directional hyperbola is 
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directed towards the driver. In this application, a 
second directional hyperbola can be directed towards the 
passenger to permit switching between driver and 
passenger during a telephone conversation as required. 

In situations in which the geometric position of the 
wanted signal source relative to the receiving 
microphones is known, geometric source separation is a 
powerful tool. The standard method of this class of "beam 
forming" algorithms is the so-called "shift and add" 
method, whereby a filter is applied to one of the 
microphone signals and the filtered signal is then added 
to the second microphone signal (see, for example, Haddad 
and Benoit, "Capabilities of a beamforming technique for 
acoustic measurements inside a moving car", The 2002 
International Congress and Exposition on Noise Control 
Engineering, Deaborn, Mi, USA, August 19-21, 2002) . 

An extension of this method relates to "adaptive 
20 beam forming" or "adaptive source separation", where the 
position of the sources in space is unknown a priori and 
has to be determined first by algorithms (WO 02/061732, 
US6, 654 , 719) . In this instance, the aim is to determine 
the position of the sources in space from the microphone 
25 signals and not, as is the case in "geometric" beam 
forming, to specify it beforehand on a fixed basis. 
Although adaptive methods have proved very useful, 
information is usually also necessary a priori in this 
case because, as a rule, an algorithm can not decide 
30 which of the detected speech sources is the wanted signal 
and which is the interference signal. The disadvantage of 
all known adaptive methods is the fact that the 
algorithms need a certain amount of time to adapt before 
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sufficient convergence exists and the source separation 
is successfully completed. Furthermore, adaptive methods 
are more susceptible to diffuse background interference 
in principle because it can significantly impair 
5 convergence. A more serious disadvantage with 
conventional "shift and add" methods is the fact that 
with two microphones, only two signal sources can be 
separated from one another and diffuse background noise 
is not attenuated to a sufficient degree as a rule. 

10 

Patent specification DE 69314514 T2 discloses a 
method of separating acoustic signals of the type 
outlined in the introductory part of claim 1. The method 
proposed in this document separates the acoustic signals 

15 in such a way that ambient noise is removed from a 
desired wanted acoustic signal and the examples of 
applications given include the speech signals of a 
vehicle passenger which can be understood but only with 
difficulty due to the general and non-localised vehicle 

20 noise. 

As a means of filtering out the speech signal, this 
prior art document proposes a technique whereby a 
complete acoustic signal is measured with the aid of two 

25 microphones, a Fourier transform is applied to each of 
the two microphone signals in order to determine its 
frequency spectrum, an angle of incidence of the 
respective signal is determined in several frequency 
bands based on the respective phase difference, which is 

30 finally followed by the actual "filtering". To this end, 
a preferred angle of incidence is determined, after which 
a filter function, namely a noise spectrum, is subtracted 
from one of the two frequency spectra, and this noise 
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spectrum is selected so that acoustic signals from the 
area around the preferred angle of incidence assigned to 
the speaker are amplified relative to the other acoustic 
signals which essentially represent background noise of 
5 the vehicle. Having been filtered in this manner, an 
inverse Fourier transform is then applied to the 
frequency spectrum which is output as a filtered acoustic 
signal . 

10 The method disclosed in DE 69314514 T2 suffers from 

the following disadvantages: 

a) The acoustic signal separation disclosed in 
this prior art document is based on completely 

15 separating an element of the originally 

measured complete acoustic signal, namely the 
element referred to as noise. In other words, 
this document works on the basis of an 
acoustic scenario in which only a single 

20 , wanted noise source exists, whose signals are, 

so to speak, embedded in interference signals 
from non-localised or less localised sources, 
in particular vehicle noise. The method 
disclosed in this prior art document therefore 

25 enables this one wanted signal exclusively to 

be filtered out by completely eliminating all 
noise signals. 

In situations where there is a single wanted 
30 acoustic signal, the method disclosed in this 

document may well produce satisfactory 
results. However, in view of its basic 
principle, it is not practical in situations 



in which not only one wanted sound source but 
several such sources contribute to the 
acoustic signal as a whole. This is the case 
in particular because, in accordance with this 
teaching, only a single so-called dominant 
angle of incidence can be processed, namely 
the angle of incidence at which the acoustic 
signal with the most energy occurs. All 
signals which arrive at the microphone from 
different angles of incidence are necessarily 
treated as noise 

Furthermore, this document itself appears to 
work on the assumption that the proposed 
filtering in the form of a subtraction of the 
noise spectrum from one of the two frequency 
spectra does not produce satisfactory results. 
Consequently, this document additionally 
proposes that yet another signal processing 
step should be performed prior to the actual 
filtering. Effectively, in all frequency 
bands, once the dominant angle of incidence 
has been determined, by means of an 
appropriate phase shift of one of the two 
acoustic signals in this frequency band to 
which a Fourier transform has been applied, 
the noise elements in the respective frequency 
band are attenuated relative to the wanted 
acoustic signals which might possibly also be 
contained in this frequency band. Accordingly, 
this document regards the filtering process 
which it discloses, in the form of a 
subtraction of the noise spectrum, as being 
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unsatisf actory in itself and actually proposes 
other signal processing steps immediately 
beforehand, which are performed by separate 
components provided specifically for this 
purpose. In particular, in addition to a 
device for subtracting the noise spectrum 
(device 24 in the single drawing appended to 
this document) , the system needs means 20 
connected upstream to effect a phase shift as 
well as means 21 to add spectra in the 
individual frequency bands after phase 
correction (see the relevant components 
illustrated in the single drawing appended to 
this document) . 

Consequently, the method and the device needed 
in order to implement it are complex. 

Accordingly, the objective of the present invention 
20 is to propose a method of separating acoustic signals 
from a plurality of sound sources and an appropriate 
device which produces output signals of a sufficient 
quality purely on the basis of the filtering step, 
without having to run a phase-corrected addition of 
25 acoustic spectra in different frequency bands in order to 
achieve a satisfactory separation, and which also not 
only enables signals from a single wanted noise source to 
be separated from all other acoustic signals but is also 
capable in principle of separately outputting acoustic 
30 signals from a plurality of sound sources without 
elimination. 

This objective is achieved by the invention on the 
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basis of a method as defined in claim 1 and a device as 
defined in claim 7. Advantageous embodiments of the 
invention are defined in the respective dependent claims. 

5 The method proposed by the invention requires no 

convergence time and is able to separate more than two 
sound sources in space using two microphones, provided 
they are spaced at a sufficient distance apart. The 
method is not very demanding in terms of memory 

10 requirements and computing power and is very stable with 
respect to diffuse interference signals. By contrast with 
the conventional beam forming process, such diffuse 
interference can be effectively attenuated. As with all 
methods involving two microphones, the spatial areas 

15 between which the process is able to differentiate are 
rotationally symmetrical with respect to the microphone 
axis, i.e. with respect to the straight line defined by 
the two microphone positions. In a section through space 
containing the axis of symmetry, the spatial area in 

20 which a sound source must be located in order to be 
considered a wanted signal corresponds to a hyperbola. 
The angle i9 0 which the apex of the hyperbola assumes 
relative to the axis of symmetry is freely selectable and 
the width of the hyperbola determined by an angle Y3db is 

25 also a freely selectable parameter. With only two 
microphones, output signals can also be created for any 
other different angles i9 0 and the separation sharpness 
between the regions decreases with the degree to which 
the corresponding hyperbolas overlap. Sound sources 

30 within a hyperbola are regarded as wanted signals and are 
attenuated with less than 3 db. Interference signals are 
eliminated depending on their angle of incidence 9 and 
an attenuation of >25db can be achieved for angles of 
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incidence 3 outside of the acceptance hyperbola. 

The method operates in the frequency range. The 
signal spectrum assigned to the one directional hyperbola 
5 is obtained by multiplying a correction function K2(xl) 
and a filter function F(f,T) by the signal spectrum 
M(f,T) of one of the microphones. The filter function is 
obtained by spectral smoothing (e.g. by diffusion) of an 
allocation function Z ( i9-i9o) and the computed angle of 

10 incidence 3 of a spectral signal component is included 
in the argument of the allocation function. This angle of 
incidence 3 is determined from the phase angle cp of the 
complex quotient of the spectra of the two microphone 
signals M2 (f , T) /Ml (f , T) , by multiplying cp by the acoustic 

15 velocity c and dividing by 2nfd, where d denotes the 
microphone distance. Having been restricted to an amount 
that is less than or equal to one on the basis of 
x=Kl(xl), the result xl=cpc/2nfd, which is also the 
argument of the correction function K2(xl), gives the 

20 cosine of the angle of incidence 3 which is contained in 
the argument of the allocation function Z(i9-i9o); in the 
above, Kl(xl) denotes another correction function. 

Brief description of the drawings: 

25 

Fig. 1 illustrates the definition of the angle of 
incidence 3 based on the positions of the two 
microphones whose signals are processed. 

30 Fig. 2 illustrates an example of an allocation 

function Z{3) with half-value width 2y3db/ which results 
in a hyperbola with the apex at 3=0. 
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Fig. 3 illustrates a hyperbola with the apex at 
i9 = i9 0/ which determines the directional characteristic of 
the source separation. Signals within the spatial area 
defined by the hyperbola are output as a wanted signal 
5 with an attenuation <3db. 

Fig, 4 illustrates the structure of the source 
separator in which the time signals of two microphones, 
ml (t) and m2(t), are transformed in a stereo-sampling and 

10 Fourier transform unit (20) to produce spectra Ml(f,T) 
and M2(f,T), where T denotes the instant at which the 
spectra occur. From the spectra, the frequency-dependent 
angle of incidence i9(f,T) as well as the corrected 
microphone spectrum M(f,T) are calculated in the &- 

15 calculating unit (30) , from which output signals (t) 
are produced in signal generators (40) for different 
directional angles i9o- 

Fig. 5 illustrates the structure of the 3 
20 calculating unit (30), in which the phase angle <p(f,T) of 
a spectral component of the complex quotient of the two 
microphone spectra Ml(f,T) and M2(f,T) is calculated, 
which then has to be multiplied by the acoustic velocity 
c and divided by 2nfd, where d notes the microphone 
25 distance. This operation gives the variable xl(f,T) which 
represents the argument of the two correction functions 
K2 and Kl . These correction functions give the corrected 
microphone spectrum M ( f , T) =M1 ( f , T) *K2 (xl ( f , T) ) and the 
variable x (f ,T) =K1 (xl (f , T) ) , from which the angle of 
30 incidence i9(f,T) is calculated by applying the inverse 
cosine function. 
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Fig. 6 illustrates a signal generator in which an 
allocation function Z(S-& 0 ) with an adjustable angle i9 0 
is smoothed by spectral diffusion to obtain a filter 
function F(f,T), which is multiplied by the corrected 
5 microphone spectrum M(f,T). This results in an output 
spectrum (f,T), from which an output signal (t) is 

obtained by applying an inverse Fourier transform, which 
contains the acoustic signals within the spatial area 
fixed by the allocation function Z and the angle 3 0 . 

10 

Fig. 7 illustrates examples of the two correction 
functions K2(xl) and Kl(xl). 

One basic principle of the invention is to allocate 

15 an angle of incidence & to each spectral component of 
the incident signal occurring at each instant T and to 
decide, solely on the basis of the calculated angle of 
incidence, whether the corresponding sound source lies 
within a desired directional hyperbola or not. In order 

20 to soften the correlation decision slightly, a "soft" 
allocation function Z ( 3 ) (Fig. 2) is used instead of a 
hard yes/no decision, which permits a continuous switch 
between desired and undesired incidental directions, 
which advantageously affects the integrity of the 

25 signals. The width of the allocation function then 
corresponds to the width of the directional hyperbola 
(Fig. 3). The complex spectra of the two microphone 
signals are divided in order to calculate, firstly, the 
phase difference cp for each frequency f at an instant T. 

30 The acoustic velocity c and the frequency f of the 
corresponding signal component are used to calculate, on 
the basis of the phase difference, a path difference 
lying between the two microphones when the signal was 
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transmitted from a point source. If the microphone 
distance d is known, the result is a simple geometric 
consideration to the effect that the quotient xl from the 
path difference and microphone distance corresponds to 
5 the cosine of the sought angle of incidence. In practice, 
due to interference such as diffuse wind noise or spatial 
echo, an assumption can rarely be made about a point 
source, for which reason xl is not usually limited to the 
anticipated value range [-1,1]. Before the angle of 

10 incidence i9 can be calculated, therefore, another 
correction factor which limits xl to said range is 
necessary. If the angle of incidence i9 (f,T) was 
determined at the instant T for every frequency f, the 
spectrum of the desired signal is obtained within a 

15 directional hyperbola with the apex at the angle 9 = 9 Q by 
a simple frequency-based multiplication by the spectrum 
of one of the microphones, in other words 
Ml (f ,T) K ( 9 (f ,T) -t9o) - Under certain circumstances, it is 
of advantage to apply spectral smoothing to K ( 9 (f , T) -i9o) 

20 before running the multiplication. Smoothing, the result 
of which is denoted by F Sq (f,T) , is obtained by applying 

a diffusion operator for example. In situations where the 
variable x used to calculate the angle of incidence lies 
outside its value range due to the effect of 

25« interference, it is of advantage to attenuate the 
corresponding spectral component of the microphone signal 
since it may be assumed that interference signals are 
superimposed. This is done by applying a correction 
function, for example, the argument of which is the 

30 variable xl . If M ( f , T) is the corrected microphone 
signal, the process of creating the desired signal 
spectrum including spectral smoothing and correction is 
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expressed by S i9o (f / T)= F^ ( f , T ) M ( f , T ) . The time signal 
(S^(t) for the corresponding directional hyperbola with 
apex angle i9 0 is obtained from S &q (t) by applying an 
inverse Fourier transform. 

5 

In other words, one basic idea of the invention is 
to distinguish noise sources, for example the driver and 
passenger in a vehicle, from one another in space and \ 
thus separate the wanted voice signal of the driver from 

10 the interference voice signal of the passenger, for 
example, making use of the fact that these two voice 
signals, in other words acoustic signals, as a rule also 
exist at different frequencies. The frequency analysis 
provided by the invention therefore firstly enables the 

15 overall acoustic signal to be split into the two 
individual acoustic signals (namely of the driver and of 
the passenger) . Then, with the aid of geometric 
considerations based on the respective frequency of each 
of the two acoustic signals and the phase difference 

20 between the output signal of microphone 1 and of 
microphone 2 associated respectively with this acoustic 
signal, it is "then only" necessary to calculate the 
direction of incidence of each of the two acoustic 
signals. Since, in a hands-free system in the vehicle, 

25 the geometry between the position of the driver, the 
position of the passenger and the position of the 
microphones is more or less known, the wanted acoustic 
signal which has to be further processed can be separated 
from the interference acoustic signal on the basis of its 

30 different angle of incidence. 

A detailed explanation of an example of an 
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embodiment of the invention will be given with reference 
to the appended drawings. 

The time signals ml(t) and m2(t) of two microphones 
5 which are disposed at a fixed distance d from one another 
are applied to an arithmetic logic unit (10) (Fig, 4), 
where they are discretized and digitised in a stereo 
sampling and Fourier transform unit (20) at a sampling 
rate f A . A Fourier transform is applied to a sequence of 

10 a sampling values of each of the respective microphone 
signals ml (t) and m2 (t) to obtain the transformed complex 
value spectrum Ml(f,T) respectively M2(f,T), in which f 
denotes the frequency of the respective signal component 
and T specifies the instant at which a spectrum occurs. 

15 In terms of the practical application, the following 
selection of parameters is suitable: f A = 11025 Hz, 
a=256, T a/2=t. If computing capacity and memory space 
permit, however, a=1024 is preferred. The microphone 
distance d should be shorter than the half wavelength of 

20 the highest frequency to be processed, which is obtained 
from the sampling frequency, i.e. d < c/4f A . For the 
parameter selection specified above, a microphone 
distance d = 20 mm is suitable. 

25 The spectra Ml(f,T) and M2(f,T) are forwarded to a 

^-calculating unit with spectrum correction (30) , which 
calculates an angle of incidence 3 (f,T) from the spectra 
Ml(f,T) and M2(f,T), which specifies the direction from 
which a signal component with a frequency f arrives at 

30 the microphones at the instant T relative to the 
microphone axis (Fig. 1). To this end, M2(f,T) and 
Ml(f,T) are subjected to a complex division. cp(f,T) 
denotes the phase angle of this quotient. In situations 
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where confusion can be ruled out, the argument (f,T) of 
the time- and frequency-dependent variables is omitted 
below. Based on the Euler formula and the arithmetic 
rules for complex numbers, the exact arithmetic rule for 
5 determining cp is as follows: 

cp=arctan ( (Rel*Im2-Iml*Re2 ) / (Rel*Re2 + Iml*Im2 ) ) , 

where Rel and Re2 denote the real parts and Iml and Im2 

10 denote the imaginary parts of Ml, respectively M2 . The 
variable xl=cpc/2nfd is obtained on the basis of the 
acoustic velocity c from the angle cp, xl also being 
dependent on frequency and time: xl=xl(f,T). In practice, 
the range of values for xl must be limited to the 

15 interval [-1,1] with the aid of a correction function 
x=Kl(xl) (Fig. 7). Taking the variable x calculated in 
this manner, an inverse cosine function is applied in 
order to calculate an angle of incidence i9 of the 
relevant signal component to be measured from the 

20 microphone axis, i.e. from the straight line defined by 
the positions of the two microphones (Fig. 1) . Taking 
account of all the dependencies, the angle of incidence 
of a signal component with frequency f at the instant T 
is therefore: i9 ( f , t ) =arccos (x ( f , T) ) . The microphone 

25 spectrum is also corrected with the aid of a second 
correction function K2(xl) (Fig. 7): 

M(f ,T) =K2 (xl)Ml (f ,T) . The purpose of this correction is 
to reduce the corresponding signal component in 
situations where the first correction function applies 

30 because it may be assumed that there is superposed 
interference which distorts the signal. The second 
correction is optional or M ( f , T) =M1 ( f , T ) may also be 
selected as an alternative; M ( f , T ) =M2 ( f , T) is also 
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possible . 

The spectrum M(f,T) together with the angle i9(f,T) 
is forwarded to one or more signal generators (40) where 
5 a signal to be output (t) is respectively obtained 

with the aid of an allocation function Z ( i9 ) (Fig. 2) and 
a selectable angle &q. This is done by multiplying every 
spectral component of the spectrum M(f,T) by the 
corresponding component of a ^o-specific filter F &o (f,T) 

10 at an instant T. (f,T) is obtained by a spectral 

smoothing of Z(&-9 0 ). This smoothing is obtained, for 
example, by spectral diffusion: 

(f,T) = Z(9 (f,T)-9 0 ) + DA 2 f Z(,9 (f,T)-,9o) . 

15 

In the above, D denotes the diffusion constant which 
is a freely selectable parameter greater than or equal to 
zero. The discrete diffusion operator A 2 f is an 
abbreviation for 

2 0 A 2 f Z{i9 (f,T)-i9„)) = (Z(i9 (f-f A /a),T)-i9 0 )-2Z< 9 <f,T)- i9 0 ) )+Z{ i9 {f+f A /a,T)-i9 0 ))/(f A /a) 2 . 

The quotient f A /a obtained from the sampling rate f A 
and number a of sampling values corresponds to the 
distance of two frequencies in the discrete spectrum. 
25 Applying the resultant filter F^ (f,T) will give a 

spectrum S i9o (f,T)= F^ ( f , T) M ( f , T) , which is transformed 

into the time signal s«9 0 (t) by inverse Fourier transform. 

The signal (t) to be output by a signal generator 

30 (40) corresponds to the acoustic signal within that area 
of space defined by the allocation function Z ( 9 ) and the 
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angle 3 0 . For the sake of simplicity , only one allocation 
function Z ( 3 ) will be used in the nomenclature selected 
for different signal generators and different signal 
generators will use only different angles 3 0 . In 
5 practice, there is nothing to say that a separate form of 
the allocation function can not be selected in each 
signal generator as well. Applying allocation functions 
permitting a decision as to different areas of space to 
which signal components belong is one of the central 

10 principles of the invention. An allocation function must 
be a direct function and appropriate functions are, for 
example, Z ( 3 ) = ( ( 1+cos 3 ) /2 ) n with a parameter n>0 . The 
spatial area in which signals are attenuated with less 
than 3db corresponds to a hyperbola with a beam angle 

15 2y3db (Fig. 3) and apex at the angle 3 0 . Accordingly, 2y3db 
corresponds to the halve-value angle of the allocation 
function Z ( 3 ) (Fig. 2), where the specified formula for 
the allocation function is y 3 db = arc cos ( 2 1_1/n -l ) . In 
these two-dimensional geometric considerations, it must 

20 be borne in mind that the actual area of the three- 
dimensional space from which acoustic signals are 
extracted with the described method is a hyperboloid of 
revolution, obtained by rotating the described hyperbola 
about the microphone axis. 

25 

Naturally, the present invention is not limited to 
use in motor vehicles and hands-free units. Other 
applications are conference telephone systems in which 
several directional hyperbola are disposed in different 
30 spatial directions in order to extract the voice signals 
of individual persons and prevent feedback or echo 
effects. The method may also be combined with a camera, 
in which case the directional hyperbola always looks in 
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the same direction as the camera so that only acoustic 
signals arriving from the image area are recorded. In 
picture-phone systems, a monitor is simultaneously 
connected to the camera, in which the microphone system 
5 can also be integrated in order to generate a directional 
hyperbola perpendicular to the monitor surface, since it 
can be expected that the speaker is located in front of 
the monitor. 

10 A totally different class of applications becomes 

possible if, instead of evaluating the signal to be 
output, the angle of incidence «9 to be determined is 
evaluated, which is then determined by averaging over 
frequencies f at an instant T, for example. This type of 

15 i9 (T) evaluation may be used for monitoring purposes if 
the position of a sound source is to be located in an 
otherwise quiet area. 

Correct "separation" of the desired area 
20 corresponding to the wanted acoustic signal to be • 
separated from a microphone spectrum need not necessarily 
be obtained by multiplying with a filter function as 
illustrated by way of example in Fig. 6, the allocation 
function of which is plotted by way of example in Figure 
25 2. Any other way of correlating the microphone spectrum 
with a filter function would be appropriate, provided 
this filter function and this correlation cause values in 
the microphone spectrum to be more intensely "attenuated" 
the farther their allocated angles of incidence 9 are 
30 from the preferred angle of incidence 9q (for example the 
direction of the driver in the vehicle) . 
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List of reference numbers 



10 Arithmetic logic unit for running the method 

steps proposed by the invention 
5 20 Stereo sampling and Fourier transform unit 

30 9 -calculating unit 

40 Signal generator 

a Number of sampling values transformed to the 

spectra Ml, respectively M2 
10 d Microphone distance 

D Diffusion constant, selectable parameters greater 

than or equal to zero 

A 2 f Diffusion operator 

f Frequency 
15 f A Sampling rate 

Kl First correction function 

K2 Second correction function 

ml(t) Time signal of the first microphone 

m2(t) Time signal of the second microphone 



20 Ml(f,T) Spectrum of the first microphone signal at the 
instant T 

M2(f,T) Spectrum of the second microphone signal at the 
instant T 

M(f,t) Spectrum of the corrected microphone signal at 
25 the instant T 

S^ (t) Time signal generated corresponding to an angle 

i9 0 of the directional hyperbola 
S^(f,T) Spectrum of the signal s5o(t) 

Y3db Angle determining the half-value width of an 

30 allocation function Z{9) 

cp Phase angle of the complex quotient M2/M1 

9 (f ,T) Angle of incidence of a signal component, 



-19- 



measured from the microphone axis 
i9o Angle of the apex of a directional hyperbola, 

parameters in Z(i9-i9 0 ) 
x, xl Intermediate variables in the 9 -calculation 
5 t Time basis of the signal sampling 

T Time basis for generating the spectrum 

Z ( 9 ) Allocation function 



